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1 Introduction 



In this work, we prove the existence and the nonlinear stability for a certain class of stationary 
solutions of the semi-relativistic Schrodinger-Poisson system in a finite volume domain with 
Dirichlet boundary conditions. Such a system describes the mean-field dynamics of semi- 
relativistic quantum particles (for instance, in the case of heated plasma), with the particles 
moving with extremely high velocities. Let us consider semi-relativistic quantum particles 
confined in a domain Q C M 3 which is an open set with a C 2 boundary and \Q\ < oo. These 
particles are interacting by means of the electrostatic field they collectively generate. In 
the mean-field limit, the density matrix describing the mixed state of the system solves the 
Hartree-von Neumann equation 

idt P (t) = [H v ,p(t)], xeQ, t>0 
—AV = n(t,x), n(t,x) = p(t,x,x), p(0) = po 

with Dirichlet boundary conditions, p(t, x, y) = if x or y G dQ, for t > 0. The single 
particle Hamiltonian is given by 

H v :=T m + V(t,x). (1.2) 

The relativistic kinetic energy operator T m := \J — A + m 2 — m is defined by means of the 
spectral calculus. In system ( 11. ip and further below, A stands for the Dirichlet Laplacian 
on L 2 (Q), and m > is the single particle mass. We refer to [3] and jl] for a derivation of 
the analogous system of equations in the non-relativistic case. Due to the fact that p(t) is a 
nonnegative, self-adjoint and trace-class operator acting on L 2 (f2), we can expand its kernel, 
for every t G R+, with respect to an orthonormal basis of L 2 (Q). We denote this kernel at 
the initial time t = by po, 

Po(x, y) = Y1 X kipk{x)^ k {y). (1.3) 
fceN 

Here {ifjk}km stands for an orthonormal basis of L 2 (Q), such that ipk\on — for all k G N, 
and the coefficients are given by 

A := {A fe } fc6N G I 1 , X k >0, J2 Xk = 1 - ( L4 ) 

fcgN 

In [1] , we showed that there exists a one-parameter family of complete orthonormal bases of 
L 2 (fl), {ipk{t)}keNi with i/)k(t)\d(i — for all k G N, and for t G M+, such that the kernel of 
the density matrix p(t), which satisfies system ( II .ip . can be expressed as 

p(t, x, y) = ^ ^k*Pk(t, x)i/j k (t, y). (1.5) 
fceN 
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As a consequence of the particular commutator structure of (11. ip (where pit) and —iHy 
satisfy the conditions of a Lax pair), the corresponding flow of pit) leaves its spectrum 
invariant. Accordingly, the coefficients A are independent of t. This isospectrality is crucial 
for the stability analysis for stationary states based on the Casimir energy method employed 
in this paper; see also [51 El [121 EES]- 

By substituting the expression (I1.5P in system (jl.ip . one can verify that the one-parameter 
family of orthonormal vectors {ipk{t)}k&N satisfies the semi-relativistic Schrodinger-Poisson 
system equivalent to (II. ip and given by 

= + V^fc, k e N, (1.6) 

at 

-AV[*]=n[% *:={Vfc}£i, (1.7) 

oo 

n[^(t,x)] = Y,Wk\ 2 , (1.8) 

k=l 

with initial data {"0fe(O)}feLi- Here, the potential function is the solution of the Poisson 
equation (11.71) . Both V[^>] and i/Jk(t), for all fcsN, satisfy the Dirichlet boundary conditions 

ij>k(t,x) , V(t, x) = 0, t > 0, Vx G dQ. (1.9) 

The global well posedness for system ( ll.6p - (ll.9l) was established in the recent work pQ. 
Analogous results were derived before in the nonrelativistic case in a finite volume domain 
with Dirichlet boundary conditions in [6], and in the whole space of M 3 in [6] and [8]. 

In this paper, we are interested in the properties of stationary states which occur when 
p(t) = f(Hy) for some function /. Substituting the latter in (jl.ip . the commutator on 
the right side of the first equation of system (11.11) vanishes, and the density matrix is time 
independent. The precise properties of the distribution function / will be discussed below. 
The solution of the Schrodinger-Poisson system corresponding to the stationary states is 

^(t ) x) = e- i "*Vfc(a?), keN, 

such that the potential function V[^>] is time independent, ^Gl are the eigenvalues of the 
Hamiltonian (}1.2j) and ipk{x) are the corresponding eigenfunctions. 

The organization of this article is as follows. In Section 2, we describe the class of 
stationary states we will study, and state our hypotheses and main results about nonlinear 
stability and existence of stationary states. In Section 3, we derive some preliminary results. 
In Section 4, we prove the nonlinear stability of the stationary states of the semi-relativistic 
Schrodinger-Poisson system via the energy-Casimir functional as a Lyapunov function (see 
the statement of Theorem 1). In Section 5 we define the dual functional and in Section 
[6] study its properties using the methods of convex analysis, and show that it admits a 
unique maximizer (see Theorem 2), which implies the existence of a stationary state for 
our Schrodinger-Poisson system. In an Appendix, we discuss generalizing the well-posedness 
result of [1] to spaces with higher regularity. 
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2 The Model and Statement of the Main Results 



The state space for the Schrodinger-Poisson system is defined as 

i 

C : = {(^, A) | ^ = {^pk}kLi C Hq (fi) n H 1 (fi) zs a complete orthonormal system in L 2 (Q), 

oo „ 

A = {A fc }£° =1 g J 1 , A fc > 0, k G N, V A fc / \Vi> k \ 2 dx < oo}, 

see p. 

In order to precisely define the class of stationary states we will study, we need to intro- 
duce the Casimir class of functions. We say that a function / : R — > R is of Casimir class C 
if and only if it possesses the following properties: 

(i) / is continuous, such that f(s) > for s < s and f(s) = when s > s , with some 
s e]o, oo], 

(ii) / is strictly decreasing on ] — oo, so], such that lim s ^_ 00 /(s) = oo, 

(iii) there exist constants e > and C > 0, such that for s > the estimate 

f(s) < C(l + s)" 5 - £ (2.1) 

holds. 

Note that So acts as an "ultra-violet" cut-off, and we can take is as large as we wish. 

Consider the quadruple (& , A , /xo, Vb) with (\P , A ) G £, /io = {/io,fc}fcLi rea ^ valued, and 
the potential function V G H^(Q) H H 2 (Q), such that the stationary Schrodinger-Poisson 
system holds 

(T m + y )^o,fe = Vo,ki>o,k, fceN, (2.2) 

oo 

- AVb = no = ^A ,A#o,fc| 2 , ( 2 - 3 ) 
fc=i 

with 

Ao, fc = /W), ken, (2.4) 

where / G C. Then, the corresponding density matrix po = f{T m + Vo) solves the stationary 
state Hartree-von Neumann equation 

[HvoiPo] = 0. 



Remark. In the nonrelativistic case, the Casimir class was defined similarly in [12] 
with the exception that the rate of decay of the distribution function / was assumed to be 



4 



smaller. A good example of / G C is the function decaying exponentially as s — > oo with the 
cut-off level So = oo. This is the Boltzmann distribution f(s) := e _/3s , (3 > 0. 

In order to prove the nonlinear stability of the stationary states, we will use the energy- 
Casimir method. This method was used in [5] for fluid problems, and in [3, [13] for studying 
stationary states of kinetic equations, in particular, Vlasov-Poisson systems. Here, we extend 
the energy- Casimir functional used in [12] to the semi-relativistic case. For / G C, let us 
define 

poo 

F(s) : = / f(a)da, s el. (2.5) 



Note that the function defined via (12.51) is decreasing, continuously different iable, nonnegative 
and is strictly convex on its support. Moreover, for s > 

F(s) < C(l + s)- 4 ~ £ . (2.6) 

Its Legendre (Fenchel) transform is given by 

F*(s):=sup AeM (A S -F(A)), s < 0. (2.7) 

We define the energy-Casimir functional for a given / as 

nc{*,\) := ^l-Afc) + W> A), A) G £. (2.8) 

k=l 

In particular, "He* is conserved along solutions of the Schrodinger-Poisson system, as a con- 
sequence of isospectrality of the flow of p(t), which is equivalent to the t-independence of A&. 
Our main results in this paper address the existence and stability of stationary states that 
are given by (12. 2D - ( 12331) . for / G C Stability is controlled by the following main theorem. 

Theorem 1. Let (ty , A ,/io, Vo) be a stationary state of the semi-relativistic Schrodinger- 
Poisson system, where 

Ao,jfe = /(/io,fc), keN 

with some f G C and (V&oiAo) G C Let {^{t),X) be a solution of the Schrodinger-Poisson 
system, such that initial datum (\l/(0), A) G C. Then, for all t > 0, the estimate 

i||n[W(t), A] - n \\%_ 1(n) < 7* c (tt(0), A) - ftc^o, A ) 

holds, such that the stationary state is nonlinearly stable. Here, H" 1 ^) is the dual of H l (Q) 
with norm ||«||^_ 1(n) = (u, (-A)" 1 ?/)^. 

To prove the existence of stationary states, we introduce the dual of the energy-Casimir 
functional. To this end, we let, for A > fixed, 

f 1 1 f c 

G(V,\,V,(t) :=^[F*(-A fc ) + A fe / [\T^ k \ 2 + V\ij k \ 2 ]dx]-- / |W| 2 dx + a[^A fe -A 

k=l J Q J Q k=l 



Here, a G R is a Lagrange multiplier. 
The dual functional to Tic is given by 

$(V,<r) :=inf* iA £($,A,l>). (2.9) 

The infimum in the formula above is taken over all A G and all complete orthonormal 
sequences ^ from L 2 (f2). Let us consider only non-negative potential functions and define 

Hl + (Q) := {V G £#(0) | V > 0}. 

The following is our main result about the existence of stationary states. 

Theorem 2. Let f G C and A > be fixed. The functional $ 

(V, a) G Hl + (Q) xl^-J f \VV\ 2 dx- Tr[F(T m + V + a)} - a A 

is continuous, strictly concave, bounded from above and — $(V, <r) coercive. There ex- 
ists a unique maximizer (Vo,o"o) of <§>(V,cr). Let {ipo,k}kLi be the orthonormal sequence 
of eigenf unctions of the Hamiltonian T m + Vo corresponding to the eigenvalues {/io,fc}£Li 
and Ao,fc := f(fJ>o,k + cr ). Then (\&o, A , [i , V ) is a stationary state of the semi-relativistic 
Schrodinger-Poisson system, where YlkLi A 0i /c = A and (\P j A ) G C 

We will prove Theorem 1 in Section 4, and Theorem 2 in Section 5. 



3 Preliminaries 

We have the following elementary lemma. 
Lemma 3. For (\&, A) E C we have 

fceN 

-^ei V^a stand for the Coulomb potential induced by n^ t \, such that 

-AV^ jX (x) = rty,\{x), ie!l; V^x) = 0, x G <91X 
T/ien V^ x e #o( fi ) n # 2 (fi). 

Proof. We easily express the norm as 



K,Alli 2 (fi) = 5^ A fc A s(l^fclM^| 2 )L2 



fc,s=l 



Here and further below, the inner product of two functions f(x),g(x) G L 2 (Q) is denoted as 
(fi9)L 2 {n) '■= Jo f(x)g(x)dx. Application of the Schwarz inequality to the right side of the 
identity above yields the upper bound 



(f>VZ 



\i>k\ 4 dx 



3^ 2 



which can be estimated from above by applying the Schwarz inequality as well. Thus, we 
obtain 

We next make use of the Sobolev inequality 

JjVf\ 2 dx>c s njf\ 6 dx) , (3.1) 



in which the constant c s is given on p. 186 of [TT] . Noting that a function compactly supported 
in the set Q can be extended by zero to the whole space of 1R 3 , we arrive at the upper bound 



\n\r ' 



c 2 



J^h \Vip k \ 2 dx < oo 



by means of the definition of the state space L given above, such that n^ t \ G L 2 (Q). Note 
that the particle density n^ t \ vanishes on the boundary of the set Q by means of formula 
(I1.8P and boundary conditions (11.9p . Therefore, AV^a G L 2 (f2). Let {//fe}feeN denote the 
eigenvalues of the Dirichlet Laplacian on L 2 (Q) and ^5 is the lowest one of them. Note that 

pi > 0, k G N. 

Since 

we have that || V^aIIl 2 ^) < - n II^aIU 2 ^) < 00 • Furthermore, V^a vanishes on the boundary 
' Mi 

of the set fi via ([EE]). ■ 

According to Theorem 1 of p], for every initial state (W(0), A) G C, there exists a unique 
strong solution of system (jl.6[) - fll.9p . where (^(t), A) G £ for all t > 0. Let us define the 
energy of a state A) G £ as 

%(*>A) :=V"A fe / |T^^ fc | 2 rfx + i / n^xV^^dx = 
k=l Jn * Jn 
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= f> fc / \ri^ k \ 2 dx + \ [ \vv iJ<x \ 2 dx, 

which is a conserved quantity along solutions of the Schrodinger-Poisson system (see Lemma 

7 of m). 

Analogously to [1] we assume that > via density arguments. To prove the nonlinear 
stability for a specified stationary state, We have the following auxiliary statements. 

Lemma 4. Let f G C. 

a) For every j3 > 1 there exists C = C(f3) G M, such that for s < we have 

F(s) >-(3s + C 

b) LetV G H^(Q) and V(x) >0/ori6fi. Then both operators f(T m + V) and F(T m + V) 
are trace class. 

Proof. The part a) of the lemma comes from the fact that function F(s) is smooth with 
the slope varying from — oo to 0, and convex; therefore, its graph is located above a tangent 
line to it. 

For the Dirichlet eigenvalues of — A on L 2 (Q), Q C M 3 , we will make use of the semiclas- 
sical lower bound 

nl > Ck 2 *, keN 

with a constant dependent on \Q\ < oo (see e.g. [10]). Since the potential function V(x) 
is nonnegative in Q by assumption, we easily estimate from below the eigenvalues \i k of the 
Hamiltonian T m + V for k G N as 

Vk> \j nl + m 2 -m> (Ck* - m)+ (3.2) 

with the right side of the inequality above positive for k large enough. For the sharp semi- 
classical bounds on the moments of Dirichlet eigenvalues to fractional powers see [H]. We 
express 

oo 

TrF(T m + V) = J2 F M<oo, 

k=l 

since F(s) is decreasing, satisfies estimate (12.61) . and the series with a general term (1 + 
(Cks — m) + )~ 4_e converges. Similarly, 

oo 

Tr/(T m + y) = ]T/(^)<oo, 

k=l 

due to the fact that f(s) decreases, obeys bound (12.1 j) and the series with a general term 
(1 + Cks — m)~ 5 ~ £ is convergent. This completes the proof of the part b) of the lemma. ■ 
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Lemma 5. Let -0 G (Q)C\H 1 (fl) with HV'lU^n) = 1, the potential function V G Hq(Q) 
and V(x) > /or Then, 

F((ifj, (T m + F)^)) < (V>, F(T m + (3.3) 

/io/ds wzi/i equality if is an eigenstate of the Hamiltonian T. m + V. 

Proof. By means of the Spectral Theorem we have 



T m + V = y^/UfcPfc, 



fc=l 

where the operators {P^}^ are the orthogonal projections onto the bound states corre- 
sponding to the eigenvalues {/Jk}kLi- Hence 

(oo 
j>*IIWllL» ( n) 
fc=i 

The right side of ( 13. 3ft can be easily written as 

oo 

fc=l 

Estimate ( j3 .3f) follows from Jensen's inequality. When if) is an eigenstate of the operator 
T m + V corresponding to an eigenvalue fik, for some k G N, both sides of ( 13. 3ft are equal 
to F(/j,k). Note that the converse of this statement does not hold in general. Indeed, let 
us consider as if) a linear combination of more than one eigenstate of the Hamiltonian with 
corresponding eigenvalues \ik located outside the support of F(s). Then both sides of (13.31) 
will be equal to zero. ■ 

The lemma below shows that a stationary solutions belong to the state space for the 
Schrodinger-Poisson system. 

Lemma 6. Let the quadruple (\l/o, A , /io, V ) satisfy equations 112.2]) , 112.3]) and \2.J$ , 
where \l/o a complete orthonormal system in L 2 (Q) and the distribution f G C. Then, 



oo „ 

y2 \k / ivv'o./oi 2 ^ < oo 

fe=l J n 



holds, such that (^ojAo) ^ £■ 

Proof. We express the following quantity using (12.21) and (12 .4p as 



00 i r 

y~] A ,fc||r^o,fc||L2(n) + / \^V \ 2 dx 



= X 0,k(( T m + V )lp 0> k, %,k)L^{n) = ^2 f(^0,k)^0,k- (3.4) 
k=l k=l 

The potential function Vq(x) > in Q since it is superharmonic by means of (12.31) . and 
vanishes on the boundary of Q. Thus, /io,& > 0, kEN and via (12.11) the right side of (13.41) 
can be bounded above by 

oo 

C ( l + ^o,fc)" 5 " e ^o,fc < oo, 

k=l 

which follows from the eigenvalue estimate (13. 2p . We have also obtained 

VV EL 2 (Q), TlVo.fc EL 2 (Q), kEN. (3.5) 

In fact, from equation ( 12. 2 p with a nonnegative potential, we easily conclude that 

i 

1 1 T£ip 0ik 1 1 2 L 2 {n) < jU 0) fe, kEN. (3.6) 

Note that the standard requirement Vq E L x (f2) (see e.g. p. 234, 245 of [11]) is satisfied here 
as well. Throughout the article the operator \p\ := v— A, which is the massless relativistic 
kinetic energy operator defined via the spectral calculus. Clearly, in the sense of quadratic 
forms, we have 

T m > \p\ - m, (3.7) 

such that (13.51) yields 

\\\p\*^o,k\\h(n) < {{T m + m)ip Qyk ,il)Q M ) L 2 {n) = m+ WT^i/jqAl^q) < °°- 
i 

Thus, ^o.jfc £ Hq (fi), kEN. Moreover, let us make use of the relativistic Sobolev inequality 
(see e.g. p. 183 of [IT]) for a function compactly supported in fl 

(f,\p\f)L>(n) >c r s \\f \\ls (n) , (3.8) 
which gives us ipo,k £ -^ 3 (^)i kEN. Hence by means of Holder's inequality, we arrive at 

jjV \ 2 \^ k \ 2 dx< (^JjV \ 6 dx^j ^J^ 0tk \ 3 dx^j <oo. 

Indeed, V (x) E L 6 (Q) due to flS3} along with (IO) . Therefore, V tp ,k E L 2 (Q), k E N. 
From equation (12. 2p we easily deduce that T m -0 O fc E L 2 (Q), k E N as well. Then the identity 

\\T m ipo,k\\l^(U) = / |V^ ,fc| 2 ^ - 27n\\TZip ,k\\h(n) ( 3 -9) 
J n 
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via (133|) yields V^o.fc G L 2 (fi) and ^o,fc G H 1 ^), k G N. By means of (12~4"1) . we have 
^o,fc > 0, k G N. Convergence of the series on the right side of (13. 4p implies 

oo oo 

^A 0ifc = ^/(/i , fc ) <oo, (3.10) 
fc=l k=l 

such that A = {Xo,k}'kLi G i 1 . Let us make use of identity (13.91) along with (13. 6p . such that 

OO n OO 

J^A , fc / |V^ ,fc| 2 ^ = ^ A ,fc{2m||T^^o,fc||i2 (n) + \\T m ipo,k\\h(n)} < 
fc=i fe=i 



< 2m^/( / u 0i fe)y[xo,fc + ^Ao i/ fc||r m ^o,fc||i2(Q). (3.11) 
fc=i fc=i 

The first term in the right side of (13.111) is finite as it was shown above. The second expression 
on the right side of the inequality above can be written via (12.21) as 

oo oo .. 

5^ \k\\no,k^o,k-V ifj ,k\\h(n) = Vfc{/^fc + 11^0^1112(^-2^ / V \tl) , k \ 2 dx}. (3.12) 
k=i k=i J Q 

Our goal is to prove that (I3.12p is convergent. Indeed, (12.11) implies 



k=i 



k=l 



due to the eigenvalue bound (I3.2p . We estimate the second term on the right side of (13.121) 
using Holder's inequality, such that 

i 

oo ( r \ 3 oo ( r ^ 

^ A o,fcllWo.Jfc||i 3 (n) < ( J n \Vo\ 6 dxj Y1 X °A J n M 3 dx 

where Vo(x) G L 6 (f2) as discussed above. Let us make use of inequalities (13. 8p and (13. 7p . 
such that 

1 1 

^2 \k\\^o,k\\h(n) < —^^oAIpI^o^^o^l^q) < — ^ A , fc {m + llTmV'o.fcll^)} < 00 



k=l 



k=l 



k=l 



due to estimates (13.101) and (13. 4p . The last term in the right side of (13.121) can be bounded 
above in the absolute value by applying the Schwarz inequality to it twice, such that 



^ Ao,fc// 0] fc||Vb^i 



0,fc||L 2 (n) 



< 



k=l 



\ k=l 
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as it was shown above. ■ 

Remark. In the stationary situation, our semi-relativistic Schrodinger-Poisson problem 
can be easily written as 

-AV Q = f(T m + V )(x,x), xeQ, 
V (x) = 0, xe dtt. 

Let us turn our attention to defining the corresponding Casimir functional for a fixed / G 
C. The following elementary lemma gives the alternative representation for the Legendre 
transform of our integrated distribution. Note that / £ C considered on the (—00, sq] semi- 
axis has an inverse 

Lemma 7. For the function F(s) defined in A2.5\) and s < we have 

F*(s) = J f-\a)da. (3.13) 

Proof. Let us define 

g(X) := Xs - F(X), AeR, s < 0, 

such that via (12. 5ft we have </(A) = s + /(A). Hence, the maximal value in the right side of 
(12. 7p is attained at A* := s) and is equal to 

POD 

cp(s) := g (X*) = f-\-s)s - / f(a)da 

with <p(0) = 0. Although / is continuous and not necessarily differentiable, we can approxi- 
mate it by a differentiable function. Let f e (s) = h f(t)dt, and 

POO 

(f€ (s) = fr 1 (-s)s- / f e (a)da. 

A direct computation using the formula for the derivative of the inverse yields (p' e {s) = 
f ( T 1 (—s). Integrating and taking the e — > limit yields <p(s) = J° s f~ 1 (o~)da. ■ 

In the next section, we prove the nonlinear stability of stationary states, using the energy- 
Casimir functional defined above. 
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4 Stability of stationary states 



In this section, we prove Theorem 1, which yields lower bound in terms of the electrostatic 
field. The auxiliary statement below is crucial for establishing this nonlinear stability result. 

Lemma 8. Let V E Hq(Q) and V > 0. (i) Then, for (\l/, A) G C, the lower bound 

ji^C-A*) + X k j[\TU k \ 2 + VW 2 ]dx j > -Tr[F(T m + V)}. (4.1) 

i 

(ii) Equality is attained for , A) = (^v,X v ), where ipv, k G Hq (Q) flif 1 (f2), fceN stands 
for the orthonormal sequence of eigenfunctions of the Hamiltonian T m + V with corresponding 
eigenvalues fiv, k and Xy, k = f(fJ>v, k ), k G N. 

Proof. According to definition (12. 7p . we have 

F*(a)>fia-F{fi), fieR, s<0, 

which easily implies 

F\-X k ) + X k fi k >-F(tJL k ), keN. (4.2) 

Now let 

fi k := / \ \rii) k \ 2 + V\i) k \ 2 \dx = (i> k , (T m + V)i/t k ), k G N, 



which proves part (i). To prove part (ii), we note that after summation, we arrive at 

+ X >< J q | \ T ^\ 2 + W fe | 2 jdx j > - J2 F ^ (T m + V)i) k )). 

Lemma 5 along with the definition of trace yields the lower bound for the right side of the 
inequality above as 

oo 

- F ( T ™ + V)ih) = "Tr(F(T m + V)). 

k=l 

Suppose that (ty, A) = (\&y, Ay), where ipv, k are eigenfunctions of the Hamiltonian T m + V 
and /ifc, which are defined above are the corresponding eigenvalues /iy; fc , k E N. Therefore, 
on the right side of the lower bound (14. ip we have 

oo 

-Tr(F(T m + V)) = -J2 F (»v, k ). 

k=l 
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Next, we use the identity Ay ifc = f(fiv, k ) = ~F\^v, k )- Then, via Lemma 7, F*'(— Xv, k ) = 
J -1 (Ay.) = fiv, k , k G N. Using the argument of Lemma 7, we arrive at 

F*(-X V J = sup AeR (— AAy fe — F(A)) = -r\\ v ^X Vt -F{r\\ Vtk )) = -\ V ^ V - F {^ k ). 

Therefore, the left side of (14. ip will be equal to — Y^k=i F{^v, k ) as wen - ■ 

Armed with the technical lemma above, we may now prove our first main statement. 

Proof of Theorem 2. Let A) G £ and the potential function V = is induced by 
this state. Then we will use the following identity for the energy of the electrostatic field 

^||n-n ||^ 1(n) = i||W-VVb||i 2(n) = i / \VV\ 2 dx+ l - f |VVb| 2 rfa;+ f V AVdx. 

By the definition of the energy-Casimir functional, this can be written as 

H c (y,\)-{jr(F*(-h) + X k f |t1^|W)-~ I \VV \ 2 dx- ( V AVdx}, 

which is equal to 

A) - { £ [F*(-A fc ) + A fc / (|T|^ fc | 2 + Uo|^| 2 )rfx] - i / |W | 2 ofe}. 



Applying first Lemma 8(i), and subsequently Lemma 8(H), we obtain that the expression 
above is bounded from above by 

n c (*, A) - { - Tr[F(T m + V )] - l - J \VV \ 2 dx} = 
= n c (y,\) - { £ [F*(-\ ,k) + A , fc / (iri^l 2 + Wo,fc| 2 )dx] - \\ \W \ 2 dx} = 

= Wc(*, A) - { j>2 [F*(-X ,k) + A , fc / \TUo,k\ 2 dx] + \ I \W,\ 2 dx} = 

= Hc(^,A)-Hc(^o,A ). 

Since the Casimir functional is constant along the solutions of the Schrodinger-Poisson sys- 
tem, which is globally well-posed (see pQ), for an initial condition (^(0), A) G £, we can use 
T-Lc(^(0), A) in the estimate above instead of ^c(^(^),A)- ■ 

After establishing the nonlinear stability of the stationary states of the semi-relativistic 
Schrodinger-Poisson system, our main goal is show the existence of such states satisfying the 
assumptions of the stability theorem. 
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5 Dual functionals 



For every distribution function f £ C we will derive a corresponding stationary state as the 
unique maxinizer of a functional defined below. Let us use the energy-Casimir functional 
from the stability result to obtain such a dual functional. Our tool below will be the saddle 
point principle. Recall that, for A > fixed 

g{V,\,V,(T) :=^[F*(-A fc ) + A fe / [\T^, k \ 2 + V\^ k \ 2 ]dx] - - / | VV\ 2 dx + a [ ^ A, - A 
fc=i J Q J Q k=i 

1 

Here as before \1/ = {ifi^kLi C Hq (f2) n H 1 ^) is a complete orthonormal system in L 2 (f2) 
and A G /+ = {(Afc) G Z 1 | A fc > 0, fe G N}. Now the function V G #o(ft) is allowed 
to vary independently of ^ and A. The parameter a G K here plays the role of Lagrange 
multipliers. The statement below demostrates how the functional defined above is related 
to our energy-Casimir functional. 

Lemma 9. For arbitrary A, a, 

00 

sup v Q(^f,X, V, a) = H C (V, A) + ° [ ^ X k - a] . (5.1) 



k=l 



The supremum is attained at V = V^. 

Proof. We express the functional defined above as 

00 L 00 

£(^,A,W) = V[F*(-A fc ) + A fc / |T r I^| 2 ^ + -A fc / \ip k \ 2 V^, x dx] + Y] A fc / V|^ fc | 2 dz- 
fe=i ^ n fe=i 

/ |V^, A | 2 rfx-i / |W| 2 dx + a[f>-A . 
Using the definition of the energy-Casimir functional (12. 8 j) we arrive at 

The expression above can be written as 

1 00 
Hc(*, X) - -||W* A - WUi^nj + a[]T A, - A 

fc=i 

which completes the proof of the lemma. 
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In the next Section, we will show that the functional $(V, o~) defined in (12.91) has a unique 
maximizer, which is a stationary state of our Schrodinger-Poisson system. Let us first prove 
the following auxiliary statement, which is the generalization of Lemma 8 above. 

Lemma 10. Let V G Hq(Q) and V > 0. Then for (\&, A) G C and a G R, the lower 
bound 

oo „ 

[F*(-\k) + \k( / [\T^ k \ 2 + V\iP k \ 2 ]dx + a)] >-Tr[F(T m + V + a)} (5.2) 
k=i ^ n 

i 

is valid. Equality in it is attained when (\&, A) = (^y, Ay), where ipv,k G i^ 2 (f^nif 1 ^), k G 
N i/te orthonormal sequence of eigenf unctions of the operator T m + V corresponding to 
eigenvalues fiv,k- Moreover, Ay,* = f(nv,k + c)> & G N. 

Proof. Let us use inequality (14. 2 j) with 

At* := / (|rf^| 2 + ^|^ fc | 2 )^ + a = (T m + V + a)ii k ), keN. 



iraVfcl 2 + V^IV*] 2 ] da; + 0") > -F((^ fc , (T m + V + a)^)), G N. (5.3) 



T m + l/ + a= / (A + o~)dE\, 



Therefore, 

F*(-A fc ) + A fc 
Clearly, 



where Ea is the spectral family associated with the Hamiltonian T m + V, such that di/ k (X) : = 
(ip k , dE\ip k ) is a probability measure for k G N. By means of Jensen's inequality 

F«Vfc, {T m + V + a)^» = F( ^°°(A + a)di/ fc (A)) < F(A + a)dz/ fc (A) = 

= (ij k ,F(T m + V + a)ij k ). 
This upper bound along with (I5.3P and summation over k G N give us the desired inequality 

i 

Then consider {"0v,fe}feLi C Hq(Q) n J/ 1 (fi) forming a complete orthonormal system in 
L 2 (f2), such that (T m + V)Vv,fc = Hv,k^v,k and Ay* = f(^v,k + cr), fc G N. In this case the 
right side of (15. 2 p is equal to 

oo oo 

- ^{F(T m + V + a)i> v ,k, M =~J2 F ^ k + a )' 
fe=i fc=i 



1(3 



We have for k G N 

F*(-X v , k ) = sup AeR (-AA y , fc - F(X)) = -r\\ v ,k)\v,k - F(f~ 1 (X Vik )), 

since it is attained at the maximal point A* := f~ l (\v,k)- The equality Xy,k = f{(J>v,k + c) 
yields f~ l (\v,k) = Hv,k + cr, such that 

F*{-\ v ,k) = -{Vv,k + cr)Xv,k ~ F(pvjt + cr)- 
A direct computation implies that the left side of (I5.2p equals to — YlT=i F{^v,k + cr)- ■ 

Armed with the auxiliary lemma above we manage to derive the expression for the dual 
functional for our problem. 

Lemma 11. The infimum in definition Ii2.9\) is attained at ^> = {ipv,k}kLi> an or ~ 
thonormal sequence of eigenf unctions of the Hamiltonian T m + V, V > corresponding to 
eigenvalues fiy,k with \y,k = /(/-ty.fc + c) f or k G N. Furthermore, the dual functional is given 
by 

$(V,a) = -- [ \W\ 2 dx-Ti[F{T m + V + 0-)} - aA. (5.4) 
2 Jn 

Proof. Let us show that the operator F(T m + V + a) is trace class. Clearly, 

oo 

Tr[F(T m + V + a)]=Y, F(fJ, V)k + a). 

k=l 

Since the potential function V > by assumption, we use inequalities (13.21) and (12.61) and 
arrive at the series with the general term (1 + Ck$ — m + a)~ A ~ e . This series is convergent. 
We conclude the proof of the lemma by referring to the result of Lemma 10 above. ■ 



6 Existence of stationary states 

In this section we prove, for each distribution function / G C and each value of A > 0, the 
existence of a unique maximizer of the functional $, which will be a stationary state of the 
semi-relativistic Schrodinger-Poisson system. 

Proof of Theorem 2. Let us first show that the bound 

Tr[F(T m + a(V 1 + a 1 ) + (l-a)(V 2 + a 2 )} < aTr[F(T m + V 1 + a 1 )] + (l-a)Tr[F(T m + V 2 + a 2 )\ 

(6.1) 

holds for any a G (0, 1) and (Vj,^) G iZj + (fi) x R, j = 1,2. Let c/> G H${Q) H H\Q) and 
II0||l 2 (Q) = 1- AVe make use of the spectral decompositions 

rco roc 

T m + V 1 = / 7 dP 7 , T m + V 2 = / pdQp, 
Jo Jo 
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where P 7 and Qp are the spectral families associated with the operators T m + V\ and T m + V 2 
respectively. This enables us to introduce the probability measures 



and write 



dv(i) := (4>,dP^(f)) L 2 {n) , dfi(P) := ((/>,dQp(l>) L 2( n) (6.2) 
F(((f>, [T m + a{Vx + a x ) + (1 - a){V 2 + <x 2 )]0) L2(n) ) = 

poo POO 

= F(a / (7 + <ri)di/(7) + (1 - a) / (/? + a 2 )d^)). 
Jo Jo 

Since F is strictly convex on its support, we obtain the upper bound for the expression above 

using Jensen's inequality as 



00 ^00 



a F( 7 + <7i)dz/(7) + (1 - a) / F(P + a 2 W(P). 
Jo Jo 

By means of definition (I6.2p we arrive at 

a(<f>, F(T m + Vt + a 1 )0) i2(n) + (1 - a)(<f>, F(T m + V 2 + a 2 )4>) L 2 {n) . 

Let {ipk}kLi be the set of eigenf unctions of the operator T m + a(Vi + G\) + (1 — a){V 2 + 02) 
forming a complete orthonormal system in L 2 (Q). Then via the argument above we obtain 

00 

F((i> k , [T m + a(Vi + <n) + (1 - a)(F 2 + a 2 )}^ k ) L 2 {n) ) < 

k=i 

00 00 
< a ^2(ip k , F(T m + Vi + ai)ifJk) L 2 (n) + (1 — a) J2^k,F(T m + V 2 + a 2 )i{j k ) 

k=l k=l 

and arrive at inequality (16.11) . Suppose equality holds. From the fact that the function F is 
strictly convex on its support we deduce that the operators T m + V\ + <J\ and T m + V 2 + a 2 
with potential functions V± and V 2 vanishing on the boundary of Q, have the same set of 
eigenvalues and the corresponding eigenfunctions are {ipk}kLi- Therefore, V\(x) = V 2 (x) in 
Q and o\ = cr 2 , and Tr[F(T m + V + a)] is strictly convex. Since — | f n \ W\ 2 dx and — crA are 
concave, we obtain that our functional given by ( 15. 4p is strictly concave. 

Then we turn our attention to the proof of its boundedness from above and coercivity. 
Obviously, by means of the Poincare inequality 

1 f i^t,|2 - ^ C l || T r||2 



with a constant C\ > 0. Let \iy be the lowest eigenvalue of the Hamiltonian T m + V. Clearly, 
we have the estimate with a trial function 6 as 



nv< / {MHr + vi^jdx, H0n i2(n) = 1. 



n 



Let us fix cf) as the ground state of the negative Dirichlet Laplacian on L 2 (Q). Then 

\\p\*i\ 2 dx = y/c;, 

where C p is the constant in the Poincare inequality. We introduce 



n 



C 2 :=\ Wdx > 0, 



which is finite. Indeed, <p G L 6 (Q) via the Sobolev inequality (I3.ip . Hence via the Schwarz 
inequality we arrive at 

r V\4>\ 2 dx<C 2 \\V\\ L 2 {n) <C 2 \\V\\ Hiin) , 



n 

such that 

Hv < ^/C~ P + C 2 \\V\\ H i {n) . 

This yields the upper bound 

m <0 < ~\\V\\^ (a) - F(VC~p + C 2 \\V\\ H . {n) + a) - aA. (6.3) 



2 

Let us use the convexity property, such that 

F(x) > -/3x + C 3 , 
where (3 > A > is large enough. This implies the inequality 

HV, a) < _^||V||^ (n) + 09 - A)a + f3C 2 \\V\\ H , {n) + p^C~ p - C s . 
A direct computation yields the estimate 

HV, a) < -Q.\\V ||^ 1(n) + C 5 + (/? - A)a + py/c 9 - c 3 . 

Let us choose j3 = 2A and define the nonnegative constant k := max{Cs + {3^fC~ v — C3, 0}. 
Hence 

m^)<-^\\Vf H}m +Ao + k. (6.4) 
Combining estimates (16.31) and (16.41) . we easily arrive at 

$(^)<-^IMlV)-AM + A;, 
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which shows that our functional $(V, a) is bounded above and — $(V, a) is coercive. There- 
fore, $(V, a) has a unique maximizer (Vo, cr ). Let the hamiltonian T m + V have the sequence 
of eigenvalues {/io,fc}fcLi an d corresponding eigenfunctions -{/0o,fc}fc=L; such that 

(T m + V Q )ip Q ,k = Ho,kipo,k, k e N 

and let Ao,fc := /(/-*o,fc + o"o)- We have 

(7) = / | VFol 2 ^ - £ / " ffA > 

such that cr = o"o is its critical point. Therefore, 

= — (Vq, o)\a=a = -A + /(A*o,* + ^0) = E A , fc - A, 

fe=l fc=l 

such that Xl^=i -V* = A. The first variation of $(V, ctq) at V = Vo vanishes as well. Thus, 
an easy computation gives us 

00 

-AV Q (x) = J2 X oM^o,k^)\ 2 - 

k=l 

By direct substitution, the functions ^(x, t) = e~ lMo ' fet ^o,fc( ;:c ) > & £ N satisfy the Schrodinger 
equation 

i^f = [T m + V }ij k , xen, t>0. 

The density matrix 

00 00 
p (t,x,y) = ^2Xo,kipk(x,t)ip k (y,t) = ^ Ao |fe -0 O|fe (x)-0o ifc (y), 
fc=i fc=i 

such that = and the particle concentration n (t,x) = po(t,x,x). 

Therefore, (^O) Ac A*0) Vo) is a stationary state of our semi-relativist ic Schrodinger-Poisson 
system. Finally, we are in position to show that (\& , Ao) e £5 which can be done analogously 
to the proof of Lemma 6 above. ■ 

We have the following result relating the functional $ and %c- 

Proposition 12. Let the assumptions of Theorem 2 hold, such that (^qj Ao> £*o> Vo) is 
the corresponding stationary state of the semi-relativistic Schrodinger-Poisson system. Then 

$(Vb,CTo)=ftc(*0,A )- 
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Proof. Note that 



$(V , a ) = -\( \VV \ 2 dx - Tr[F(T m + V + a )\ - a A 
1 Jn 



and 

oo 



^c(^o,A ) = ^F*(-A , fc ) + ^A , fc / \Ti^ k \ 2 dx+- \VV \ 2 dx. 

1 1 7 1 J Q f2 



fc=l fc = l 

By means of Lemma 10 



oo r i 

[F*(-Ao, fc ) + Ao, fc ( / [IT^I 2 + V \il>o,k\ 2 ]dx + ct )] = —Tr [F(T m + V + a )], 
fc=i ^ n 

which implies the statement of the proposition. ■ 
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Appendix: Higher regularity 

In this Appendix, we extend the global well-posedness result of [T] to spaces with higher 
regularity. For seM, let 

i 

C s : = {(*, A) | \& = {ipk}kLi C Ho (^)n^(n) is a complete orthonormal system in L 2 (fi), 

oo „ 

A = {A fe }~ =1 G l\ X k > 0, k G N, VaJ \(-A) s / 2 x/; k \ 2 dx < oo}. 

k=i Jn 

Note that in the state spaces defined above the Dirichlet Laplacian was used and no additional 
boundary conditions were required. We introduce inner product (•, -)u%(n) which induces the 
generalized inhomogenous Sobolev norm 

oo 

\Mni(n) ■= EA fc ||0 fe ||| s(n) )i, 

k=l 

and define the corresponding Hilbert space 

-Hl(Q) := {$ = {<M£°=i I fa e H*(tt) n H S (Q), VfceN, ||$||« i( n) < oo}. 
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We also introduce the generalized homogenous Sobolev norm 

oc 
k=l 

Clearly, H^H^^) ~ ll^ll^(Q)- Furthermore, it follows from the Poincare inequality that 
ll^llw^o) ^ This implies the following result. 

Lemma A.l. ForT G the norms ||r||^s(Q) and ||r||^ s ^ are equivalent. 

We know from pQ that, for $ G 

ll^]|ki(n)<ll^ll^/ 2(n) ll*lki(Q) 

and that 

\\F[V] - F[*]\\«i<p) $ (11*11*1(0) + H $ %n))H* " nuiw- 
We have the following inequalities in spaces of higher regularity. 

Lemma A.2. Let V, $ G s > 2. Then 

\\Fm\Hi(U) < \\n 2 nl -^ Q) ll*ll«i(n) 

and t/iat 

- F[$]|| Hi(n) < (||¥|& J(n) + ||*|ft. (n) )||* - *||« i(0 ,. 



We start by proving the first inequality. The second inequality follows using a similar 
analysis. 

WFmWnttn) < \\mfn m = E A fc A,(V[*]^, (-A) S F[*]^) 

fc,«>0 

= (-i) s E ^ A fc AK^[^]^,9 2a y[vi/]a 2s - 2Q ^) 

o<H<sfc,z>o 

0<|«|<s fc>0 

s( E E Afe ii aQy ^ii^wii 5S "^ii^(«) + E Afe ii F ^ii^(«)ii 5S ^ii^(«)) 2 

0<|a|<s-l k>0 k>0 

< ( E ll^[*] IU6 ( n) ||*|| . s _ a+ i + ll^ll^^l^lk^Q)) 2 

0<|a|<s-l 
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where we have used the generalized Leibnitz rule on the third line, Holder's inequality on 
the fourth and fifth lines, and the Sobolev inequality 



6 



< 



HP (ft) 



on the sixth line. It follows from the Sobolev inequality that 

wvmwUu) < iibr 1/2 ^[*]ii| 2( n)- 

Furthermore 

12 I 



<\\niy 2f j\\pr 1/2 n[n\^a), 



where we have used Holder's inequality in the first line, and the Sobolev inequality on the 
second line. This yields |||p| -1/ ' 2 n[W]|| i 2(n) < ||\I/|| 2 . 1/2 , and hence 

\\v[m^a)<\\n 2 ^ ay 

We now estimate ||<9 a V A (-0)|| L 6(Q), < \a\ < s — 1. 



0<\P\<\a\-l k,l>0 

L ' , fc|lL6(f7)ll^'^llL3 (f7) 



< E E A ^H 5a "^lli e (n)ll^" : 



0<\P\<\a\-l k,l>0 

E||\T/||2 ll\]>ll 2 

0<|/3|<|a|-l A 



Combining the above inequalities yields 



\m\\ 2 n m $ ( E E ll*H^(n)ll*IU^ ro JI*L^^ + ll*ll^,J*lk(n)) S 



1^(0)11-11^^11-11^^ f ll-ll^/ 2(n) ll-imi(n)) 2 

0<|a|<s-l 0<|/3|<|<*|-1 — — 



< ll^ll 4 , ll^ll 2 

and hence 

ll^[*]||«i(n)< ll*ll«.-i(n)ll*lki(n). 
To prove the second inequality, note that 

\\F[v\ - Fim-Hiia) = \\vm(y - $) + - ^[$])*ii« i(n) 

< - $)n w . (n) + \\(vm - ^])$n« s(n) . 
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An analysis similar to the proof of the first inequality yields the local Lipschitz continuity 
and completes the proof of the lemma. ■ 

Lemma A. 2 together with the fact that the operator T m generates the group e~~ tTmt , t G R, 
of unitary operators implies local well-posedness in C s . Furthermore, we know from [1] that, 
for \&(x,0) G \\W(x,t)\\ n i(a) is bounded for all times. 

Theorem A. 3. For every initial state (^(x, 0), A) G C s , s > 2, there is a unique mild 
solution ty(x,t), t G [0, oo), of U.6]) - /[l7^) with (ty(x,t),X) G C s , which is also a unique 
strong global solution in C\(Q). 

The proof follows from the blow-up alternative and the first inequality in Lemma A. 2. 
The mild solution of the Schrodinger-Poisson system (ll.6|) - (ll.8p . given by 

qr(t) = e~ iTmt ty{0) + e~ lTmt [ e iTmt ' F[^f(t')]dt' , (6.5) 

Jo 

which implies 

ll*Wlk(n)<||*(0)|| W + / \\FMt')}\\ Him dt'. 

-Jo 



If ~ 1) ^ follows from Lemma A. 2 that 

\Mt)\\ Wxm <c 1 + c 2 [ \\y(t')\\ Wxm dt'. 

Jo - 

By Gronwall's lemma, 

IW)lk(^) <C ie C2 *, t G [0, T'] C [0,T). 

Since \\^{t) \\n\(n) * s bounded for all times, it follows by induction on s that ||\l/(t)||-^s(o) is 
bounded for all times. ■ 
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